A Multi-Parameter Investigation of Gravitational Slip 
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A detailed analysis of gravitational slip, a new post-general relativity cosmological parameter char- 
acterizing the degree of departure of the laws of gravitation from general relativity on cosmological 
scales, is presented. This phenomenological approach assumes that cosmic acceleration is due to 
new gravitational effects; the amount of spacetime curvature produced per unit mass is changed in 
such a way that a Universe containing only matter and radiation begins to accelerate as if under the 
influence of a cosmological constant. Changes in the law of gravitation are further manifest in the 
behavior of the inhomogeneous gravitational field, as reflected in the cosmic microwave background, 
weak lensing, and evolution of large-scale structure. The new parameter, ujo, is naively expected 
to be of order unity. However, a multiparameter analysis, allowing for variation of all the standard 
cosmological parameters, finds that ujo = 0.09^q '59 (2c) where ujo = corresponds to a ACDM 
universe under general relativity. Future probes of the cosmic microwave background (Planck) and 
large-scale structure (Euclid) may improve the limits by a factor of four. 



I. INTRODUCTION 

Cosmic acceleration 0, E[ can be caused by new fluids, 
new theories of gravity, or some admixture of both 
This uncertainty places a premium on descriptions of the 
so-called "dark physics" which remain useful across dif- 
ferent models and in spite of varying assumptions. In the 
case of new fluids (dark energy), the literature chooses 
to speak in terms of the equation of state w and its 
derivative In the case of new gravitational physics, 
the model-independent lingua franca is the relationship 
between the Newtonian (ip) and longitudinal (<fi) grav- 
itational potentials. The potentials, implicitly defined 
through the perturbed Robertson- Walker metric 



ds 2 = a 2 [-(l + 2vb)dT 2 + (1 - 2(j))dx 2 



(1) 



are most familiar for their roles in Newton's equation, 
x = —Vip, and the Poisson equation, V 2 = 4TrGa 2 8p 
under general relativity (GR). 

The gravitational potentials are equal in the presence 
of non-relativistic stress-energy under GR. Alternate the- 
ories of gravity make no such guarantee. Scalar-tensor 
d Q and f(R) theories @, i, [|, b rancworld scenarios 
such as Dvali-Gabadadze-Porrati gravity [IE EH EH , and 
massive gravity [IE 03] all predict a systematic difference 
or "slip" , so that <f> ^ ip in the presence of non-relativistic 
stress-energy. Efforts to develop a paramctrized-post- 
Fricdmannian (PPF) framework to phcnomenologically 
describe this behavior are .just as prolific: Refs. EE EE 
El EE EE EE EH El HI Hi all offer parametrizations 
quantifying the departure from <f> = ip due to new gravi- 
tational effects. We choose to work with the parametriza- 



tion proposed in Ref. [if 

vb = 
vo(z) = 



[l + w{z)]4> 
wq{1 + z)- 3 . 



(2) 
(3) 
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We assume the existence of a theory of gravitation that 
leads to an expansion history that is indistinguishable 
from that produced by a spatially-flat, ACDM scenario 
with density parameters Q m and Q\ = 1 — O m . This 
assumption is not essential, but it allows our analy- 
sis to focus solely on PPF effects. Our naive expec- 
tation is that 137 ~ Q\/£l m by today. [Note that we 
have changed our notation, having previously defined 

w{z) = zu (n A /n m ){i + z)- 3 .} 

The departure from GR kicks in only when the 
cosmic expansion begins to accelerate. Daniel et al. 
[2^ | (hereafter DCCM) discuss the compatibility with 
other parametrizations (especially that of Ref. [23|) and 
compare the implications of wq ^ to data from 
the Wilkinson Microwave Anisotropy Probe (WMAP) 
EH , the Canada- France-Hawaii Telescope Le gac y Su rvey 
(CFHTLS) 0, and various galaxy surveys 0, EE Ej- 
We expand upon their analysis in this work by perform- 
ing a full likelihood analysis of the cosmological parame- 
ter space. 

The previous work by DCCM considered the effects of 
modified gravity on cosmological perturbations in a one- 
parameter context: i.e., "how does the the new (modi- 
fied gravity) parameter affect cosmological data when all 
other parameters are held fixed (at the WMAP 3 year 
maximum likelihood values)?" They used a modified 
version of the Boltzmann code CMBfast [3(3| to evalu- 
ate the effect of wq on the cosmic microwave background 
(CMB) anisotropy, matter power spectrum, weak lensing 
convergence correlation function, and galaxy-CMB cross- 
correlation power spectrum. While this analysis was use- 
ful for testing for the existence of PPF effects, the results 
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glossed over degeneracies that exist between zuq and tra- 
ditional cosmological parameters. DCCM's Figure 9 al- 
ready demonstrates a potential degeneracy between zuo 
and ag. Identifying further degeneracies and more rigor- 
ously motivating the possiblity of non-zero zuq requires 
analysis across the full cosmological parameter space. 

In the following, we present the results of a likelihood 
analysis based on a Monte Carlo Markov chain sampling 
of the space of cosmological parameters. The parameters, 
{Ctbh , n c h 2 , 9, r r ;, n Sl A s , A^Zi t^oli ar e respectively the 
baryon density, cold dark matter density, the ratio of 
the sound horizon to the angular diameter distance, the 
optical depth to last scattering, the scalar spectral index, 
the amplitude of the primordial curvature perturbations, 
and a normalization parameter for the SZ effect. These 
are the standard parameters in the convention used by 
the publicly-available code CosmoMC [3l| . 

We gene rate our Markov Chains using CosmoMC 
[32l l33l J3~i] with modules added to calculate likelihoods 
based on the weak lensing [H|,[36[ and galaxy-CMB cross- 
correlation spectra [37||. The CMB data and likelihood 
code comes from the WMAP team's 5-year release [381 ]. 
Supernova data comes from the Union data set produced 
by the Supernova Cosmology Project [3{|. The weak 
lensing data comes from the CFHTLS weak lensing sur- 
vey 

To help understand our results, we present a closed sys- 
tem of ordinary differential equations describing the evo- 
lution of </> and the matter overdensity S under zuq 7^ 0. 
These results imply a correction to the Poisson equation 
that was neglected in DCCM. Section |TT] presents these 
equations and uses them to describe the dependence of 
the large-angle CMB anisotropy on wq. Section [TTT] dis- 
cusses the modifications made to the public CosmoMC 
codes to implement Eq. @. Section HVl presents the like- 
lihood contours found from our Markov chains. Section 
Ivl makes an attempt at forecasting the results of future 
experiments. We conclude in Section [VTl 



II. EVOLUTION OF PERTURBATIONS 

The procedure for evolving the matter and metric per- 
turbations is as follows. We assume that the perturbed 
stress-energy tensors for all matter and radiation are con- 
served independently of the theory of gravitation: 



5e-06 



0. 



(4) 



We next impose the relationship given by Eq. ([2]) be- 
tween potentials <fi and ip, which upon translation into 
synchronous gauge implies an evolution equation for the 
metric variable a = (h + 6r))/2fc 2 : 

a = -{2 + m)Ha + {l + m)ji-12irGa 2 (p + p)(j/k 2 . (5) 

Here, a dot indicates the derivative with respect to con- 
formal time, h and r\ are the synchronous-gauge metric 
perturbations, TL = a/a is the conformal-time Hubble pa- 
rameter, and a is the shear in a fluid with mean density p 
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FIG. 1: The potentials cj> and tj) are shown versus the scale 
factor, for different values of zuo. The blue, dark curves are cj>, 
whereas the green, light curves are tp. Note that they behave 
oppositely; when <j) becomes shallower, tf) becomes deeper, and 



vice versa. 



and pressure p. (We use the same notation as Ref. |4l|.) 
We further assume that there is no preferred reference 
frame introduced by the new gravitational effects; there 
is no "dark fluid" momentum flux or velocity relative 
to the dark matter and baryon cosmic rest frame. This 
condition is imposed by enforcing the same perturbed 
time-space equation as in GR, 



k 2 i) = 4irGa 2 (p + p)6, 



(G) 



where 6 is the divergence of the velocity field in a 
fluid with mean density p and pressure p. Satisfying 
this equation automatically means that Bertschinger's 
consistency condition, that long-wavelength curvature 
perturbations should evolve like separate Robertston- 
Walker spacetimes, is satisfied [lj|. The model of zu(z) 
plus the three Ens. (BE1 ) close the system of equa- 
tions. (See Refs. [lg, [22| for further details.) In order 
to study the late-time behavior of the system of equa- 
tions, we may neglect the shear and velocity perturba- 
tions and express the evolution equations in conformal- 
Newtonian/longitudinal gauge as 

4> = -(3 + zv)H<p-mH(t)- {l + m)(H 2 +2H)<f>, (7) 



, k 



+ (l+m)H4> 

1-n/n 2 



(8) 



where S is the matter density contrast. 

Consider the behavior of an overdense region, 5 > 0, as 
it evolves from early times when GR is valid to late times 
when new gravitational effects characterized by w be- 
come important. At early times, when the Poisson equa- 
tion is valid, < for the overdensity. While the expan- 
sion is matter-dominated, the potential remains static. 
However, at late times, with the onset of cosmic acceler- 
ation, the potential begins to evolve. In the case of GR, 
4> > so the potential is stretched shallower. The density 
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FIG. 2: The matter density contrast is shown versus the scale 
factor, for different values of wq. The time evolution is ob- 
tained by integrating Eqs. Q and with initial conditions 
4> = -10~ 5 , 4> = 0.0 for k = O.OlMpc -1 . Positive (negative) 
values of wo enhance (slow) the growth of density perturba- 
tions. 
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FIG. 4: The CMB quadrupole moment is shown versus wo- 
As explained in the text, the quadratic dependence can be 
understood in terms of the influence of wo on the ISW ef- 
fect. (Reproduced from DCCM with our new normalization 
Eq. ©.) 
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FIG. 3: The degree of deviation from the Poisson equation 
versus scale factor is shown for different values of wo ■ Because 
w is scale-independent, so too is the ratio —k 2 cj)/{^-KGa 2 8p). 
For positive (negative) wo, a given <j> corresponds to a larger 
(smaller) density contrast than in GR. 



contrast 5 continues to grow via gravitational instability, 
although the rate of growth is slowed. The evolution of 
(j) can be understood in terms of a competition between 
the expansion diluting the matter density and stretching 
<f> shallower, and the accretion of matter sourcing and 
deepening <fi. In GR, the accelerated expansion upsets the 
balance in favor of dilution, so that 4> becomes shallower 
and 6 grows more slowly. When wo ^ 0, the compe- 
tition between effects changes. Numerically integrating 
Eqs. ([7]) and ([8]), we find that wq > causes 4> to become 
even shallower, yet the density contrast grows faster, as 
illustrated in Figs. [T] and O This seems counter-intuitive, 
since the shallower potential should provide weaker at- 



traction for the accretion of surrounding matter. In the 
case wq < 0, the potential <fi becomes more negative 
or deeper, and the density contrast grows more slowly. 
Likewise, the deeper potential should provide greater at- 
traction. But here the difference between <fi and ip is 
important. As seen in Fig.[TJ the potential (f> grows shal- 
lower (deeper) for vjq > (< 0). However, the potential 
responsible for acceleration ip = (1 + zu)^ behaves oppo- 
sitely, becoming deeper (shallower) . Hence, the competi- 
tion swings in favor of increased clustering over dilution 
by the expansion. 

The new behavior of 4> and 5 implies a correction to the 
Poisson equation. As seen in Fig. [3J for zuq > (< 0), 
the density contrast grows more (less) rapidly and the 
potential <fi becomes shallower (deeper) , so that the ratio 

T = -k 2 <j>/(4irGa 2 6p) (9) 

grows smaller (larger) . This suggests that we can restore 
the Poisson equation by introducing a time-dependent 
gravitational constant G e // = GT, whence —k 2 (j> = 
AirG e f fa 2 8p. Note that G e // is not a free function, but is 
determined by Eqs. (J^HHl) ■ Because we have chosen vj to 
be scale- independent, G e // is too. A different strategy, 
whereby the time- and space-dependence of G e f / is im- 
posed separately [2(|, will not necessarily satisfy Eqs. (gj- 

We can use this new understanding to explain the cu- 
rious behavior of the large-angular scale CMB anisotropy 
spectrum. The effect of wo ^ on the low I moments 
of the CMB anisotropy is not monotonic, as seen in 
Fig. 0J The cause is the suppression of the integrated 
Sachs- Wolfe effect (ISW) at vjq ~ 1. If the gravitational 
potentials in the Universe are evolving with time, CMB 
photons will lose less (more) energy climbing out of po- 
tential wells than they gained falling in, resulting in a net 
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FIG. 5: The conformal time derivative of the gravitational 
potential <f> is shown versus scale factor, for different values of 
wo. Initial conditions are the same as in Fig. [2] The potential 
well is decaying when ^/\4>iHo\ > 0, and is deepening when 
negative. 



blue (red) shift as the potentials shrink (grow). This is 
the ISW effect whereby time-evolving gravitational po- 
tentials contribute to the moments of the photon distri- 
bution function, 6;(fc, rf), via 



dr(<t>(k, r) + r))ji (fc(r - r)) exp[-r ri (*)] . (10) 



(See equation 8.55 of Ref. [42|.) Here ji is a spherical 
Bessel function of the first kind, r is the conformal time, 
To is the conformal time at z = 0, and r r i(z) is the optical 
depth to rcdshift z. The strength of the ISW effect is 
determined by the sum (jj + ip, which, using Eqs. (OOi is 
given by 



= 0(2 + vj) + 3(f>Hzu. 



(11) 



Again consider the evolution of an overdensity 6 > 
with cj> < 0. In GR, the sum is positive, <j> + ip > 0. 
When mo < 0, the second term in Eq. (fTTj) is always 
positive. The first term is generally subdominant, since 
\<p\ < H<j), as can be inferred from Fig. [5] Therefore 
nj < enhances the ISW effect. When wq > 0, there is a 
competition between the first and second terms; the first 
term is positive, whereas the second term is negative, The 
first term always wins, but at some intermediate value 
of zuq the two terms nearly cancel, thereby suppressing 
the ISW effect relative to the case with vuo = 0. This 
explains the dip in the quadrupole moment versus wq, 
seen in Fig. 0J 



III. IMPLEMENTATION 

The modifications of the Monte Carlo Markov chain 
software CosmoMC to allow for wq ^ proceed al- 
most identically to the modifications made to CMBfast 



for DCCM, with a few differences. To compare the pre- 
dictions of our model with weak lensing data we adapt 
the weak lensing module provided by Refs. [H, [36|. We 
modify it to assess the likelihood in terms of the vari- 
ance of the aperture mass (Eq. 5 of [26[) with a full 
covariance matrix [43| . Because we probe weak- lensing 
at non-linear scales, we calculate the power spectrum of 
the lensing potential by extrapolating the linear matter 
power spectrum P$ to non-linear scales and using the 
relationship @ between the matter overdensity 5 and 
the gravitational potential <p to find the non-linear P^. 
Whereas CMBfast calculates the non-linear matter power 
spectrum from the phenomenological fit of Peacock and 
Dodds [3], CosmoMC (having been built around the 
code CAMB [33[) uses Smith et aVs fit [45| (see their 
Appendix C). Smith et al. express their fit as a non- 
trivial function of the linear power spectrum and Cl m . 
This function assumes the ACDM relationship between 
f2 m and perturbation growth. Gravitational slip alters 
this relationship, as discussed above in section[TTl There- 
fore, to adapt Smith et al.'s fit to the case wq ^ 0, we 
use the phenomenological relationship (DCCM equation 
24) 

^m|ro =0 = ^m|ro #0 + 0.13tI7o— — (12) 

"A 

to find a zuq = 0, ACDM model with a similar growth his- 
tory to our -cuq 7^ model and use that value of TO | roo= o 
in Smith et al.'s equations (C18). Eq. (|12|) breaks down 
for O m | roo ^o ^ 0.15, but this region of parameter space 
is excluded to at least 2a (see Figure [6]). A second Cos- 
moMC run with a more accurate fitting function yielded 
identical results to those obtained using Eq. (fT2"|) . 

This is not a precise method for determining the non- 
linear power spectrum in the presence of gravitational 
slip. Precision would require examination of N-body sim- 
ulations which, unfortunately, implies assumptions about 
what alternative theory of gravity we are constraining. 
Recently, much work has been done attempting to calcu- 
late the non-linear power spectrum directly, without the 
aid of an N-body simulation. Crocce and Scoccimarro 
propose to expand the non-linear power spectrum as a 
Taylor-like sum 



(13) 



where the different orders of P$ are derived from a di- 
agrammatic scheme similar to Feynman diagrams [46| . 
They find that the resulting sum (flU)) is much better be- 
haved than results derived from perturbation theory (see 
their Figure 1). Matarrese and Pietroni [47j use the for- 
malism of renormalization group theory to derive a gen- 
crating functional for the different orders of P$. Taruya 
and Hiramatsu adapt methods from the statistical stud- 
ies of fluid instabilities to separate out and solve for the 
cross- mode interactions in 5 [48| . All of these methods 
yield better agreement with the results of N-body simu- 
lations than standard perturbation theory in the case of 
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w = (see Figure 2 of Ref. [4j|, Figure 8 of Ref. (471, and 
Figure 3 of Ref. [H(|). Work has already begun adapt- 
ing them to alternative gravity theories. In Ref. [511 ] - 
Koyama, Taruya, and Hiramatsu extend the method of 
Ref. [48[ to include f(R) and DGP gravity theories by 
assuming that they can be approximated with a Brans- 
Dicke scalar-tensor theory on sub-horizon scales. Hira- 
matsu and Taruya (Hoj also try to encompass modified 
gravity theories by parametrizing them in terms of their 
implied effective Newton's constant G c g = TG (see equa- 
tion [H] of the present work) . Following their lead, it 
should be possible to adapt the non-linear power spec- 
trum calculations of Ref. [H| - or even [46[ and [13] - to 
account for model-independent gravitational slip. Such 
a calculation is beyond the scope of this work. Given 
the relatively well-behaved regions of parameter space 
allowed by experiments (see Section IIVI below) , we do 
not expect this limitation to significantly influence our 
findings. 

To incorporate the galaxy-CMB cross-correlation, we 
use the module written by Ho et al. [37| ■ Modifications 
for tz7 ^ enter as modifications to the cj> + if) power 
spectrum (see section II of (37|). 
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(14) 

Note that equation (27) of DCCM neglected the factor F, 
defined in Eq. ([9]), to correct the Poisson equation. The 
corrected weak lensing statistics show the same qualita- 
tive behavior as in Figure 10 of DCCM. However, large 
values |n7o| 3> 1 have a weaker effect on the amplitude of 
the convergence spectrum. 



IV. RESULTS 

The results of our multiparameter investigation are 
shown in Figs. [3 Fig. [6] shows the 68% and 95% 
contours in (fl m ,zuo) space marginalized over all other 
parameters. Fig. [7] shows the same contours in (<7s,zdq) 
space. Red (smaller) likelihood contours were generated 
using all available data sets (WMAP 5 year [38|], Super- 
nova Union 13911 , CFHTLS [26| , and the galaxy surveys 
selected by |33|). Blue (larger) contours were generated 
using only the WMAP 5 year data. For each set of con- 
straints, we generated four independent Markov Chains. 
We achieved convergence by running the calculations un- 
til the statistic 1 1 — R\ was much less than unity, where R 
is Gelman and Rubin's potential scale reduction factor, 
defined as the ratio of the variance across all of the chains 
to the mean of the variance of each individual chain eval- 
uated for the least converged parameter. [3l|, [52|, [53| . 
Our conclusions are three-fold: 

• Present cosmological data constrains gravity to 
agree with GR, assuming a background evolution 
consistent with ACDM. 
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FIG. 6: The 68% and 95% likelihood contours in the zno — fim 
parameter space are shown. The blue contours are based on 
CMB data alone. The red contours add weak lensing, type 
la supernovae, and galaxy-CMB cross-correlation data. 
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FIG. 7: The 68% and 95% likelihood contours in the vjo — era 
parameter space are shown. Shading is the same as in Fig. [6] 
Note that the addition of large-scale structure data breaks the 
degeneracy in zuo — a$. 



• Very negative values of zuq are ruled out. This 
should not be surprising, since a sign difference 
between the longitudinal and Newtonian gravita- 
tional potentials would mean that test particles are 
repelled by overdense regions. 

• Large-scale structure data (in our case, weak lens- 
ing and the galaxy-CMB correlation) are critical to 
constraining zuq. 

The effects described in section iHl mean that any CMB 
anisotropy spectrum can be reasonably well approxi- 
mated (modulo a normalization) by two possible values 
of mo. DCCM Figure 1 showed that tzjq ^ has no 
effect on the shape of higher I multipolcs within linear 
theory. This explains the double-peaked likelihood curve 
in DCCM Fig. 3 and the broad blue contours in Figs. [S] 
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FIG. 8: The projected 68% and 95% likelihood contours in 
the oto — O m parameter space are shown. The yellow contours 
are based on mock Planck data. The green contours add mock 
weak lensing data. The underlying model is assumed to be 
vdq = with f2 m = 0.26. The current constraints are shown 
for reference. 



and [7] in this work. Fortunately, the effect of wq 7^ 
on cosmic structure is monotonic in the range of interest 
(as discussed in DCCM), so that only one value of wq is 
maximally likely for any given realization of weak lens- 
ing and galaxy-CMB cross-correlation data, hence the 
smaller red contours in Fig. [5] and [7] Marginalizing over 
all other parameters, the WMAP 5 year data alone gives 
1.7±2.o(2(t). Including supcrnovac, weak lensing, 
and the galaxy-CMB cross-correlation data improves the 
constraint to wo = O.CW^Q'gg (2a). Table [J presents the 
marginalized la limits on the other cosmological param- 
eters of note. 



V. FORECASTS 

It is useful to ask how much better our constraints will 
be under future experiments. We generate two mock data 
sets - one simulating the results of the upcoming Planck 
CMB experiment, the other simulating the results of a 
future weak lensing survey, modeled after the proposed 
ESA experiment Euclid - and feed them into our modi- 
fied CosmoMC. 

To simulate Planck data we use a fiducial model given 
by the best fit parameters of WMAP [38| with noise prop- 
erties consistent with a combination of Planck 100-143- 
217 GHz channels of HFI f54|; in this case we fit also for 



B-modes produced by lensing of the CMB (see Ref. 55 1) 
and we use the full-sky likelihood function given in 56 1 . 

To simulate weak lensing data, we generate a mock 
convergence power spectrum P K (l) (equation (2) of 
Ref. [26|]) corrected for alternative gravity as in Eq. (fT4")) . 
We generate data in bins of size A; = 1 for 2 < I < 100 
and A; = 40 for 100 <l < 2980. We simulate the (la) 



FIG. 9: The projected 68% and 95% likelihood contours in 
the wq — as parameter space are shown. Shading is the same 
as in Fig. [H] 



errors as (Eq. 11 of Ref. [57| ) 



07 



^/(2/(2Z + l))/(A,/ Bky )(P s (Z) + <7 £ 7%al), 



taking a t = 0.25, n ga i = 35(arc minute) and / s k y = 
0.48, consistent with values projected for ESA's Eu- 
clid experiment (Table 1 of Ref. These assump- 
tions will give us a tighter constraint than if we had 
used SNAP/JDEM parameters, since SNAP/JDEM has 
a smaller / s k y by a factor of 10 (59|. We fit the rcdshift 
distribution of sources n(z) from a mock data set based 
on Eq. 14 of Ref. [26| with parameter values taken from 
their Table 1. The la errors in our mock n(z) are reduced 
from actual values [H| by a factor of 1/V2. The likeli- 
hood relative to the mock weak lensing data is calculated 
as a simple \ 2 (i- e -i we assume that the covariance matrix 
is diagonal) . This is a safe assumption according to [60l | . 
Figs. [5] and O show the resulting likelihood contours. 

Looking at the Planck-only (yellow) contours, we see 
the weakness of using CMB measurements alone to con- 
strain vdq, as a bimodal distribution is obtained once 
again. We also see more clearly in Fig. [9] the degen- 
eracy between wq and ag as normalization parameters 
(one can interpret the effect of vjq on 5 in Fig. [5] as a 
renormalization of the matter power spectrum). Since 
weak lensing statistics depend sensitively on the power 
spectrum normalization, they once again break the de- 
generacy. Marginalizing over all other parameters, the 



mock datascts give the constraint w$ — —0.07 



+0.13 
-0.16 



(2*), 



a factor of ~ 4 improvement over the current constraint. 



VI. CONCLUSIONS 

If we are justified in describing the background evolu- 
tion by a ACDM universe, then the results illustrated in 
Figs. [5] and [7] do not appear to indicate a significant de- 
parture from GR. In fact, these results conflict with our 
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parameter 




UJQ = 


WMAP 5-year 


£l b h 2 


n n99K9+ 00059 

u - uzzuz -0. 00058 


02264+ ' 00058 

U.UZZO*_ 000 5 7 


0.02273 ± 0.00062 




0.1167 ±0.0026 


0.1170 ±0.0016 


0.1109 ±0.0062 


e s 


1 n417+ 0029 


l-0419±g;8Sg 


1.0400 ± 0.0029 


Tti 


0.085 ± 0.016 


0.087l° ;S; 


0.087 ±0.017 


n s 


0.964 ± 0.014 


0.965 ± 0.014 


0.96312:811 


Qa 


0.712 ± 0.014 


n 7i n+ 012 

U. 1 lU_o.011 


0.742 ± 0.030 


as 


0.842 ± 0.014 


0.844 ±0.015 


0.796 ± 0.036 


h 


0.696 ± 0.014 


0.695 ±0.013 


n 71 Q+O-026 



TABLE I: Marginalized (lcr) constraints for cosmological parameters resulting from Monte Carlo Markov chain analysis. The 
left and center columns are generated using all available data sets (CMB, weak lensing, supernovae, and galaxy-CMB cross- 
correlation). The left column is generated allows vjq, to vary. The center column fixes zuo = 0. Because our constraint on zuq is 
consistent with ujq = 0, we find little difference between the two columns. The right column shows the constraints reported by 
the WMAP team in Ref. [Hf based on just the WMAP 5-year data. The principal improvements from adding supernova, weak 
lensing, and galaxy-CMB cross-correlation data lie in constraining (a result of adding the supernovae) and as (a result of 
adding weak lensing). 



naive expectation that wq ~ Q\/Q m . However, these re- 
sults allow the ratio of <j> to ip to vary by order unity from 
the predictions of GR at the present epoch. (Weaker con- 
straints yet result if the redshift dependence of w(z) is 
allowed to vary; see Ref. [55[.) These are not very tight 
constraints. As shown in Sec. [Vj it seems likely that ex- 
periments already under consideration will give us much 
tighter constraints on parametrized-post-Friedmannian 
departures from GR in the near future. If, indeed, fu- 
ture constraints improve, we may need to reconsider the 
assumption of homogeneous zu. 

Throughout this paper wc neglect any possible scale- 
dependence of n7. This simplifying assumption seems jus- 
tified given the absence of any significant departure from 
GR. Were we to see evidence of a departure from GR, the 
onus would be on us to demonstrate the new theory's con- 
sistency with solar system-scale tests, all of which prefer 
GR to one part in 10 5 (e.g. Ref. [6l|). Beyond this exper- 
imental evidence, we expect that w should be scale de- 
pendent simply due to the differing evolution histories of 
sub- and super- horizon perturbation modes. Other work 
has already attempted tackling this expectation. Hu and 
Sawicki implement a scale-dependent gravitational slip, 
based on the behavior seen in f(R) models of gravity. 
[l§ | . Afshordi et al. offer a scale-dependent parametriza- 
tion of —(ill — (j))/(cj) + ip) designed to be consistent with 
higher-dimensional generalizations of DGP gravity [62| . 
Though they find that their parametrization is capable of 
describing effects qualitatively consistent with tensions in 
current data sets, none of those tensions is strong enough 
to warrant a detection of alternative gravity. Koivisto 



and Mota explore a different set of new gravitational ef- 
fects by supposing that dark energy is an imperfect (non- 
zero shear) fluid (r33j . Shear a, like gravitational slip, 
affects the space-space, off-diagonal perturbed Einstein 
equation, k 2 (<p — ip) = \2nGa 2 p(l + w)a. The imperfect 
fluid introduces a dark flow, however, so that the gravi- 
tational effects arc not fully equivalent to the results of 
gravitational slip. Like the present work, they find that 
data cannot yet definitively rule in or out the interesting 
regions of their parameter space. Specifically, they find 
that the effect of non-zero shear on the CMB anisotropy 
spectrum is weaker than the effect of w demonstrated in 
DCCM. 

We have also shown that the modification of the Pois- 
son equation follows uniquely from the assumptions of 
our model: the enforced relationship between (p and ip, 
stress-energy conservation, and the absence of a preferred 
frame indicated by a "dark flow" . This must be taken 
into account when conducting future tests of GR on cos- 
mological scales. 
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